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Comment on "Phase Diagram of a Disordered Bo- 
son Hubbard Model in Two Dimensions" 

In a recent Letter p| (see also Q) the authors pre- 
sented numerical evidence supporting an idea of a direct 
transition between the superfluid (SF) and Mott insulat- 
ing (MI) phases in the disordered Bosonic system, and 
even studied non-trivial properties of the multicritical 
line where SF, MI and the Bose Glass (BG) phases meet. 
The results were obtained from Monte Carlo simulations 
of the (2+l)-dimensional classical loop-current model |3J 
with the lattice action 


, vj=o 

rr 


.J 2 (v,t)-2^ + v{v))J t (v,t) 


(1) 


where r, r are spatial and imaginary time coordinates, 
and J(r, r) are integer current vectors with zero diver- 
gence. The spatial disorder potential v(r) is uniformly 
distributed on the interval (—A, A). 

Here we prove that all of the above mentioned con- 
clusions are incorrect and originate from doing simula- 
tions for too small system sizes (the maximum system 
size considered in [l| was L x x L y x L T = 14 x 14 x 20), 
and ignoring the rigorous theorem of 0, Q saying that 
for disorder A larger than the half-width of the energy 
gap in the ideal Mott insulator (we denote it as E g ) the 
system state is compressible, i.e. it is BG. Indeed, in the 
infinite system one can always find arbitrary large regions 
with the chemical potential being nearly homogeneously 
shifted downwards or upwards by A. There is no energy 
gap then for the particle transfer between such regions, 
and they can be doped with particles/holes. [The con- 
jecture is that the MI-BG transition is exactly at the 
upper boundary of the theorem, A = E g (K).] We note, 
that the theorem is based on rare statistical fluctuations, 
and for A — > the distance between regions contribut- 
ing to non-zero but exponentially small compressibility is 
diverging exponentially. 

In Fig. 1 we plot our data for E g (K) |6j along with 
the critical values of disorder for the superfluid-insulator 
transition. Since A is always larger than E g (K c ), the 
transition is always from SF to the compressible insu- 
lating phase, or BG. Accordingly, the multicritical line 
where SF, MI, and BG phases meet does not exists for 
non-zero A. 

The diverging distance between rare statistical real- 
izations of disorder is the major problem in interpreting 


Monte Carlo data for small system sizes. Recently devel- 
oped Worm algorithms |6j, Ijj allow simulations of system 
sizes as large as 160 x 160 x 1000, but for A = 0.4 even this 
is barely enough to resolve very small but finite compress- 
ibility at the critical point 6] . Clearly, for smaller system 
sizes and smaller values of disorder the short length-scale 
behavior will mimic a direct SF-MI transition observed 
in Refs. 0,0- Similar arguments were used previously 
in the criticism of the direct SF-MI transition in one di- 
mensional systems 0. 



FIG. 1: The phase diagram of the disordered, commensurate 
loop-current model Q. All error bars are of order 10 -3 and 
smaller than point sizes; data points for the SF-BG line were 
taken from: A = 0.4 Q, Q, H, A = 0.2 A = 0.1, 0.5 [this 
work] . 
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